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ESTIMATION OF SNOW TEMPERATURE AND 
MEAN CRYSTAL RADIUS FROM REMOTE 
MULTISPECTRAL PASSIVE MICROWAVE MEASUREMENTS 


A. T. C. Chang 

Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION 

A statistical analysis technique has been developed in this paper to estimate the snow 
temperature and mean crystal radius of snowfields over Greenland and Antarctica. Recent 
experimental work (Reference 1) has demonstrated that variations in crystal size and the 
physical temperature of snow field observations from space give large variations in the mi- 
crowave brightness temperature. Since the observed brightness temperature for several 
wavelengths show distinct responses to snow parameters, the Scanning Multichannel Micro- 
wave Radiometer (SMMR) measurements on board a spacecraft would appear to be 
applicable in extracting the snow temperature and mean crystal radius profiles. The SMMR 
experiment, to be launched on board the Nimbus-G and Seasat-A spacecraft, will make 
observations in wavelengths of 0.8, 1 .4, 1 .7, 2.8, and 4.6 cm. The estimated error for the 
parameters derived from this technique is approximately 1.5 K for temperature and 0.001 
cm for mean crystal radius in the presence of 1 K rms noise for each SMMR channel. 

MICROWAVE EMISSION FROM SNOWFIELDS 

The Mie scattering theory (Reference 2) was used to explain the phenomena of electromag- 
netic wave propagation within a snowfield. Since snow and ice fields generally consist of 
closely packed non-spherical particles, two assumptions had to be made before applying the 
Mie scattering theory. It was first assumed that the shape of snow and ice particles were 
spherical. Secondly, it was assumed that the particle scatters incoherently and was indepen- 
dent of the distance between scatterers. Under these two assumptions, the extinction and 
scattering cross section for a snowfield could be calculated by utilizing the Mie scattering 
theory. These quantities were subsequently used in solving the radiative transfer equation 
numerically. 

In the microwave region, the Rayleigh-Jeans approximation to the plane radiative law applies 
quite well. Therefore, the intensity of thermal radiation is directly proportional to the 



brightness temperature. The radiative transfer equation (Reference 3) for a scattering media 
may be written as 


dT B (0,Z) j 

COS9 ~~~dZ Tabs (Z) (T B (0,Z) - T (Z)) 

r ( 1 ) 

= T B (0 S , Z) F ( 0 , 0 S ) sin 0 S d 0 S -y sca (Z) T B (0,Z) 

2 J 0 


where 0 is the propagative direction, T B (0, Z) is the brightness temperature at depth Z in 
the 0 direction, T(Z) is the thermodynamic temperature of the absorbing media, F(0,0 S ) is 
the scattering phase function, 7 a b s (Z) and 7 sca (Z) are the volume absorption and scattering 
coefficient at depth Z, respectively. 

In Camp Century, Greenland, the approximate crystal size for the top 20 m firn cover was 
0.5 mm (Reference 4) which is approximately one-twentieth the size of the shortest SMMR 
wavelength. The size of the crystalis well within the range of the Rayleigh equation so that 
the Rayleigh scattering calculations can be used instead of the Mie scattering equation to 
economize on computer time. The radiative transfer equation is then solved by the invariant 
imbedding technique (References 5 and 6). 

When the snow depth is great, the reflection and emission from the ground surface, below 
the snowfield, can be neglected and the brightness temperature can be expressed as 

C°° - J z y(Z) sec 0dZ 

T b (0,O)= T(Z) 7 (Z)e J ° dZ. (2) 

Jo 


The factor multiplying T(Z) in the integrand of equation (2) is by definition the temperature 
weighting function. It represents the effective radiative depth within the medium of the 
snowfield. Figure 1 shows the calculated temperature weighting function for each of the 
SMMR channels. 


STATISTICAL RETRIEVAL TECHNIQUE 

Snow temperature, at discrete levels, can be estimated from SMMR measurements by corre- 
lation with the temperature over the weighting function layer. Since the correlation is not 
exact, the derived temperature at discrete levels from SMMR data will not be as accurate as 
the average temperature over the weighting function layer. 

The method used for obtaining the temperature and mean radius profiles from SMMR mea- 
surements are essentially a regression analysis of snowfield temperature and mean crystal 
radius, and the numerically-computed snow emission. In a statistical sense, it was attempted 
to minimize the mean-square deviation of the estimated profile from the a priori profile. 
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The statistical retrieval method used in this study is similar to those employed by C. D. 
Rodgers (Reference 7) and N. E. Gaut (Reference 8). Let us consider an ensemble consis- 
ting of sets of snow temperature profiles T(Z) or mean crystal radii r(Z) whose elements 
give the temperature at discrete depths of the snowfield and the corresponding ten channels 
of brightness temperature 0(T B ). It is assumed that the snow temperature T (Z) and mean 
crystal radius r(Z) are related to the brightness temperature 0 (T B ) by matrix D and D' 
respectively 


T (Z) 

= D 

' 0(T b ) 

(3) 

r(Z) 

= D' 

■0(T„). 

(3') 


The measured data vector 0(T B ) used for determining the snow temperature is defined as 
the measured brightness temperature minus the ensemble mean 


1 

T Biv ' T B 1V 
T B ih - T B 1H 


0 (T B ) = 


(4) 


T B 


5V 


T B 


t b 5H - t b 


5V 


5H 


where Tg n is the brightness temperature from SMMR channel n and an overbar indicates its 
ensemble mean. The inferred snow temperature and mean crystal radius vector were chosen 
to have the elements corresponding to the depth levels of: (1) surface, (2) 1 meter, (3) 4 
meters, (4) 7 meters, and (5)10 meters of the snowfields which approximately relate to the 
peak of the weighting functions. 


Matrix D was chosen to minimize the expected error between the SMMR determined profile 
and the actual temperature profile, given by 


C(T,0) C 1 (0,0) 

(5) 

C(r,0) C _1 (0,0) 

(S') 
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where C is a correlation matrix, the ij th elements is (X,Y) = XjYj and the Cy can be 
calculated from snow statistics. When C (0,0) is singular or nearly singular its inverse is un- 
defined. Under this condition, the matrix inverse calculated for this matrix C' 1 (0,0) will 
introduce considerable errors. The singularity problem was avoided by including a random 
noise in the simulated brightness temperatures. This method of handling singularities is sim- 
ilar to J. W. Waters et al. (Reference 9). To derive the D or D' matrix from equation (5) or 
(5'), it is necessary to have the C(T,0) and C(r,0) in cross-correlation matrices of the snow 
temperature and mean crystal size with the measured data vector, and C' 1 (0,0) the inverse 
of the auto correlation matrix of the data vector. Since it does not appear to be practical to 
collect the necessary ensemble by measurements, the two vectors were obtained from a priori 
measured field data and the radiative transfer calculations. The snow data used were com- 
piled by H. J. Zwally (Reference 10) with measurements taken from several stations in 
Greenland and Antarctica. In Table 1 the mean annual surface temperature (T ) and the 
scattering coefficient per unit length (7 ) are listed for each station location. The physical 
temperature profile of the snowfield can be expressed as a function of T , Z, and t (time of 
year). 


T(Z) = T m - 15e-0-3Z cos [0.99(t-84)-(97+20Z)] (6) 

The maximum surface temperature occurs at t = 0 and the minimum surface temperature at 
t = 365/2. By utilizing equation (6) and the scattering coefficient, it is possible to determine 
the expected brightness temperature emerging from the snowfield for various snow conditions. 

Table 1 


Location and Mean Annual Surface Temperature and Scattering Per Unit Length ( Zwally. 1 977) 


Location 

T m *(K) 

7s( m_1 ) 

South Pole, Antarctica 
90°S 

222 

0.222 + 0.00863Z 

Plateau, Antarctica 
79°15'S; 40°30'E 

216 

0.220+ 0.0275Z 

Camp Century, Greenland 
77°1 l'N; 61°10'W 

249 

0.1 63 + 0.0647Z 

Byrd, Antarctica 
79°59'S; 120°0l’W 

245 

0.152 + 0.0968Z 

Inge Lehmann, Greenland 
77°57'N; 39°1 l'W 

243 

0.162 + 0.1 178Z 

Site 2, Greenland 
76°59'N; 56°04'W 

249 

0.0921 + 0.021 2Z 

South Ice, Antarctica 
81°57'S;28°50'W 

242 

0.0422 + 0.0805Z 


*Mean annual surface temperature. 
Z is snow depth in meters. 
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To demonstrate the behavior of such an estimation scheme, an ensemble of eighty -four dif- 
ferent snow temperature and mean crystal radius profiles were used to generate the combi- 
nation of snow crystal and temperature profiles from the surface, down to 25 m in depth, 
in 1 m intervals. In Figures 2 and 3, the ensemble mean profile is shown along with the 
standard deviation for the ensemble at each level. Also plotted is the residual standard de- 
viation, where the rms noise is 1 K for each SMMR channel. Figures 2 and 3 demonstrate 
that the variance is reduced, in this case, where the SMMR data is used in retrieval as com- 
pared with the a priori statistical ensemble. 

To test the behavior of the retrieval scheme, Figures 4 and 5 show the results of retrieval by 
using statistical data. For the Camp Century, Greenland site (77° 1 l'N, 61°10'W), both the 
mean crystal profile and snow temperature profile are similar to the ensemble mean. Under 
these conditions, the retrieved parameters fit very well with the known parameters. In the 
South Pole site (90° S), the snow crystal radius profile differs considerably from the ensemble 
mean. Hence, the retrieved crystal radius profile does not fit as well as in the previous case. 
For the temperature profile, the retrieved data still gives good estimates as the variability of 
the snow temperature is much less than the mean snow crystal size. 

CONCLUSION 

A statistical retrieval technique has been developed to estimate the snow temperature and 
mean crystal radius profiles from SMMR data. This technique is used in the study of the 
snow temperature and mean crystal radius profiles of snowfields over Greenland and Ant- 
arctica where the snow truth is difficult to obtain by conventional methods. 

The accuracy of retrieval is influenced by the variance of the parameters in the chosen sta- 
tistical ensemble, and the success of this method depends on having a good estimate for the 
snow temperature and crystal size correlations with the brightness temperatures for the lo- 
cation in question. The ensemble used consists of data from stations with different snow 
accumulation characteristics at various times of the year. The detailed structure of the snow- 
field, such as the ice layers, are not included as further study is required to determine from 
the geographical locations which method might be reliably used. In addition, the effect of 
the atmosphere on the retrieved data should be taken into account. 

The computational method presented in this paper may also provide a basis for utilizing 
SMMR data to further develop a new retrieval technique for eventual incorporation into 
satellite weather and climate data collecting systems. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland January 1978 
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Figure 4. Retrieval data from Camp Century, Greenland (summer profile). 
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